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Fermionic bright soliton in a boson-fermion mixture 
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We use a time-dependent dynamical mean-field-hydrodynamic model to study the formation of 
fermionic bright solitons in a trapped degenerate Fermi gas mixed with a Bose-Einstein condensate 
in a quasi-one-dimensional cigar-shaped geometry. Due to a strong Pauli-blocking repulsion among 
spin-polarized fermions at short distances there cannot be bright fermionic solitons in the case of 
repulsive boson-fermion interactions. However, we demonstrate that stable bright fermionic solitons 
can be formed for a sufficiently attractive boson-fermion interaction in a boson-fermion mixture. We 
also consider the formation of fermionic solitons in the presence of a periodic axial optical-lattice 
potential. These solitons can be formed and studied in the laboratory with present technology. 

PACS numbers: 03.75.Lm, 03.75.Ss 



I. INTRODUCTION 

There cannot be an effective evaporative cooling lead- 
ing to a trapped quantum degenerate Fermi gas (DFG) 
due to a strong repulsive Pauli-blocking interaction at low 
temperature among spin-polarized fermions p|. How- 
ever, it has been possible to achieve a DFG by sym- 
pathetic cooling in the presence of a second boson or 
fermion component. Recently, there have been success- 
ful observation 0, 0, OJ H and associated experimental 
H H and theoretical H ©, GjJ GH E2 

studies of mix- 
tures of a trapped DFG and a Bose-Einstein condensate 
(BEC) by different experimental groups 10, HUH in the 
following systems: 6 ' 7 Li 0, 23 Na- 6 Li @ and 87 Rb- 40 K 
0, |ij . Also, there have been studies of mixtures of two- 
component trapped DFGs in 40 K p| and 6 Li Q atoms. 

The formation and collapse of a DFG in a boson- 
fermion mixture 87 Rb- 40 K have been observed and stud- 
ied by Modugno et al. [HElE] Although, the fermion- 
fcrmion interaction at short distances is repulsive due to 
strong Pauli blocking and hence incapable of leading to 
collapse, a sufficiently attractive boson-fermion interac- 
tion could overcome the Pauli repulsion and could result 
in a collapse of a DGF. Bright solitons in a BEC are 
formed due to an attractive nonlinear atomic interaction. 
As the interaction in a pure DFG at short distances is 
repulsive, there cannot be bright solitons in a DFG. 

In this paper we study the possibility of the forma- 
tion of stable fermionic bright solitons in a mixture of 
a DFG with a BEC in the presence of a sufficiently at- 
tractive boson-fermion interaction which can overcome 
the Pauli repulsion among fermions. The formation of 
a fermionic soliton is related to the fact that the sys- 
tem can lower its energy by forming high density regions 



(the solitons) when the attraction between the bosons 
and fermions is large enough to overcome the Pauli re- 
pulsion in the DFG and any possible repulsion in the 
BEC. In particular we consider the formation of fermionic 
bright solitons, which can freely move in the axial direc- 
tion, in such a mixture for a quasi-one-dimensional cigar- 
shaped geometry using a coupled time-dependent mean- 
ficld-hydrodynamic model where the bosonic component 
is treated by the mean-field Gross-Pitaevskii equation 
[T^l and the fermionic component is treated by a hy- 
drodynamic model This time-dependent mean- 

field-hydrodynamic model was suggested recently by the 
present author to study the collapse dynamics of 
fermions and is a time-dependent extension of a time- 
independent model used for the stationary states by Ca- 
puzzi et al. 

Bright solitons are really eigenfunction of the one- 
dimensional nonlinear Schrodinger equation. However, 
the experimental realization of bright solitons in trapped 
attractive cigar-shaped BECs has been possible under 
strong transverse binding which, in the case of weak or 
no axial binding, simulates the ideal one-dimensional sit- 
uation for the formation of bright solitons. The dimen- 
sionless nonlinear Schrodinger (NLS) equation in the at 
tractive or self-focusing case 
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with four parameters. The parameter B represents the 
amplitude as well as pulse width, v represents veloc- 
ity, the parameters 5 and a are phase constants. The 
bright soliton profile is easily recognized for v = 5 = as 
\u(x, t)\ = v2~B sech|%vi?1 ■ There have been experimen- 
tal [lj| and theoretical [lfj studies of the formation of 
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bright solitons in a BEC. In view of this, here we study 
for the first time the possibility of the formation of a sta- 
ble fermionic bright soliton in a boson-fermion mixture. 

In recent times there have been routine experimental 
studies on the formation of BEC in the presence of a 
periodic axial optical-lattice potential These leads 
to a different condition of trapping from the harmonic 
trap and generates a BEC of distinct modulation. Hence 
we also consider in this paper the modulations of the 
fermionic bright solitons in the presence of an optical- 
lattice potential. 

In Sec. II we present the time-dependent mean-field 
model consisting of a set of coupled partial differential 
equations involving the bosonic and fermionic probabil- 
ity densities. In the case of a cigar-shaped geometry 
with stronger radial trapping, the above model is reduced 
to an effective one-dimensional form appropriate for the 
study of bright solitons. In Sec. Ill we present our results 
for stationary axially-free fermionic bright solitons as well 
as those formed on an axial periodic optical-lattice poten- 
tial. We also demonstrate the stability of the bright soli- 
tons after a perturbation is applied. The bright solitons 
are found to execute stable breathing oscillation upon 
perturbation. Finally, a summary of our findings is given 
in Sec. IV. 



II. NONLINEAR 
MEAN-FIELD-HYDRODYNAMIC MODEL 

The time-dependent Bose-Einstein condensate wave 
function ^(r, t) at position r and time t may be described 
by the following mean-field nonlinear Gross-Pitaevskii 
equation p| 



2m B 



Vs(r) + gBBriB 



¥ fl (r,f)= 0,(2.1) 



with normalization J dr\9g(r, t)\ 2 = N B . Here tub is the 
mass and Nb the number of bosonic atoms in the con- 
densate, tib = |^s(r, t)\ 2 is the boson probability den- 
sity, qbb = ^h 2 a B B/rnB the strength of inter-atomic 
interaction, with asB the boson-boson scattering length. 
The trap potential with axial symmetry may be written 
as Vg(r) = |msw 2 (p 2 + v 2 z 2 ) where u> and vuj are the 
angular frequencies in the radial (p) and axial (z) direc- 
tions with v the anisotropy parameter. The probability 
density of an isolated DFG in the Thomas- Fermi approx- 
imation is given by [To| 



n F 



[max(0,{e F -V F (r)})] 3 / 2 
A 3 / 2 



(2.2) 



where A = ?i 2 (67r 2 ) 2 / 3 /(2m F ), e F is the Fermi energy, 
m F is the fermionic mass, and the function max denotes 
the larger of the arguments. The confining trap poten- 
tial Vp(r) — ^rapuo F (p 2 + v 2 z 2 ) has axial symmetry as 
the bosonic potential Ve(r), where u F is the radial fre- 
quency. The anisotropy parameter v will be taken to be 



zero for axially-free solitons in the following. The num- 
ber of fermionic atoms Np is given by the normalization 
/ drn F (r) = N F . 

We developed a set of practical time-dependent mean- 
field-hydrodynamic equations for the interacting boson- 
fermion mixture starting from the following Lagrangian 
density [13 
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(2.3) 



where g BF = 2nh 2 aB F /inR with the boson-fermion re- 
duced mass wir = rriBmp I '(niB + m F ) , where asp is the 
boson-fermion scattering length. 

The terms in the first round bracket on the right-hand 
side of Eq. I|2.3[) are the standard Gross-Pitaevskii terms 
for the bosons and are related to a Schrodinger-like equa- 
tion [Tflj . However, terms in the second round bracket are 
derived from the hydrodynamic equation of motion of the 
fermions [Tl| . Hence, the second kinetic energy term has 
a different mass factor 6m F and not the conventional 
Schrodinger mass factor 2ms as in the first integral. Fi- 
nally, the last term in this equation corresponds to an 
interaction between bosons and fermions. The interac- 
tion between fermions in spin polarized state is highly 
suppressed due to Pauli blocking and has been neglected 
in Eq. (|2.3|) and will be ne glec ted throughout this paper. 

Recently, Jezek et al. [l8j used the Thomas-Fermi- 
Weizsacker kinetic energy term Tp of fermions in their 
formulation which, in our notation, will correspond to 
a fermionic kinetic energy of H 2 \\7 r ^/n F \ 2 / (9m F ) in Eq. 
(j2.3|l in place of the present term Ti 2 \ SI ry /np\ 2 / (Qm F ). 
This kinetic energy term contributes little to this prob- 
lem compared to the dominating 3^4|nF| 5 ^ 3 /5 term in Eq. 
(|2.3|) and is usually neglected in the Thomas-Fermi ap- 
proximation. However, its inclusion leads to an analytic 
solution for the probability density everywhere [18j. For 
a discussion of these two fermionic kinetic energy terms 
we refer the reader to Rcfs. |lllll8l 

With the Lagrangian density (|2.3|) , the Euler-Lagrange 
equations of motion become [l2j : 



ih 



d h 2 vl 



dt 2m b 



V B {r) + gBBriB 



gs F n F 



*a(r,t) = 0, 



(2.4) 



ft 2 V 2 
6m F 



V F (r) + A\n F 



2/3 



3 



9BF n B 



^n F (r,t) = 0. 



(2.5) 



When the nonlinearity in Eq. (|2.5[1 is very large, the 
kinetic energy term in this equation can be neglected and 
the time-independent stationary form of this equation 
becomes 



up 



[max(0, {e F - V F (r) - cjbf n B })} 3/2 

A 3/2 



(2.6) 



which is the generalization of Eq. (|2.2I) in the presence 
of boson-fermion coupling. Equation i|2.6|) has been used 
by Modugno et al. |l(J for an analysis of a DFG-BEC 
mixture. In actual experimental condition the nonlinear- 
ity in Eq. (|2.5|) is quite large and Eq. (|2.6(l is a good 
approximation. 

For the study of bright solitons we shall reduce Eqs. 
(|2.4() and (|2.5() to the minimal one-dimensional form un- 
der the action of stronger radial trapping. The one- 
dimensional form is appropriate for studying bright soli- 
tons in the so-called cigar-shaped quasi-one-dimensional 
geometry where v < < 1. For radially-bound and axially- 
free solitons we eventually set v = 0. In this case 
the dynamical equations can be reduced to strict one- 
dimensional coupled NLS equations without any trap. 
We perform this reduction below where we take Vb (r) = 
Vf(v) = |r«BW 2 (p 2 + v 2 z 2 ) which corresponds to a re- 
duction of ojf and uu>f in Vf(t) by a factor y ms/mi? 
as in the study by Modugno et al. |lfl | and Jezek et al. 

m 

For v = 0, Eqs. (|2.4() and (|2.5(l can be reduced to an 
effective one-dimensional form by considering solutions 
of the type V B (r,t) = 4> B (z,t)t(} < ^' (p) and ^fn F (r,t) = 
4> F (z,t)^p\p), where 



nh 



exp 



M^p 7 



(2.7) 



and i = B,F represents bosons and fermions and Mb = 
m B and M F = ^/Smsmp. The expression (|2.7|) corre- 
sponds to the respective ground state wave function in 
the absence of nonlinear interactions and satisfies 



2m,B 



6m F 



(0) , 1 2 2 ,(0) 



(2.8) 
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+ -jm B u p ipp' 



m B % ,(o) 
om F 



(2.9) 



with normalization 2-k J °° ^^{p)^ pdp = 1. Now the dy- 
namics is carried by <f>i{z,t) and the radial dependence 
is frozen in the ground state tp\ (p) . The factorization 
of ^>b and yJriF above follows from the structure of the 
mathematical equations 1)2. 4f> and (|2.5|l . Although np 
gives the probability density of DFG it may not be to the 
point to associate ipF and 0f to physical fermionic one- 
particle wave functions. The true fermionic wave func- 
tion has the form of a many-particle Slater determinant. 



Nevertheless, the functions 4> F and (f>p could be regarded 
as mathematical functions related to fermionic density 
y|. In the quasi-one-dimensional cigar-shaped geometry 
the linear fermionic and bosonic probability densities are 
given by \4>f(z, t)\ 2 and |</>b(z, t)\ 2 , respectively. 

Averaging over the radial mode ^P{p), i.e., multiply- 
ing Eqs. I|2.4|) and (|2.5|l by ipf 3 (p) and integrating over 
p, we obtain the following one-dimensional dynamical 
equations |2rH : 



in— 



dt 2ttib dz 2 
+ Fbb\4>b\ 2 + Fbf\4'f\ 



a h 2 d 2 

dt 6m f dz 2 



<t> B (z,t)=0, (2.10) 
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1 4/3 



F, 



BF\<PB\ 



4> F {z,t) = Q, (2.11) 



where 
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Fbf = Qb F 
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pdp Mbfuj 

= 9bf ~j , (2.13) 
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(2.14) 



In Eqs. (|2.1U|) and (|2.11|) the normalization is given by 
/_ \4>i( z it)\ 2 dz = iVj. In these equations we have set 
the anisotropy parameter v — to remove the axial trap 
and thus to generate axially-free quasi-one-dimcnsional 
solitons. 

For calculational purpose it is convenient to reduce 
the set <|2.1(Jll and (|2.11|) to dimensionless form by in- 
troducing convenient dimensionless variables. Although 
the algebra is quite straightforward, the expressions be- 
come messy with different factors of masses. As we shall 
not be interested in a particular boson-fermion system in 
this paper, but will be concerned with the formation of 
fermionic bright solitons in general, we take in the rest 
of this paper mj = imp = m( 87 Rb), where m( 87 Rb) 
is the mass of 87 Rb, and whence m B = 3mp / '4, Mb = 
Mp = ms, and Mbf = fn B /2. In the two experimental 
situations of Refs. 0,0 m B ~ 3m f. 

In Eqs. (|2.1U|) and H2.11fl . we consider the dimension- 
less variables r = tu/2, y — z/l, \i = \J (jTNijfiii with 
I = yjh/ (wtob), so that 

d_ 

'8t 



d 2 
dy 2 



i— _ — + Nbb \^ b \ 



Nbf \Xf\ 



Xb(v,t) = 0, 



(2.15) 



4 



d_ 

"97 



~ t + Nfb \xb\ 



N F f \xf 



dy 2 

4/3 



Xf(v,t) = 0, 



(2.16) 



where N BB = 4clbbNb/1, N B f = 8a B FN F /l, N F b = 
8a B F N B /l , and N FF = 9(6irN F ) 2/3 /5. In Eqs. (£Ej) 
and l|2.16l) . the normalization condition is given by 
f-oo \Xi(y, T )\ 2 dy = 1. Equations 1|2.15[) and (|2.16|) are 
the coupled one-dimensional NLS equations describing 
the formation of solitons in the DFG-BEC mixture in a 
cigar-shaped quasi-one-dimensional geometry. 

In Eqs. (OS) and (|^TB|) the term N F f\xf\ 4/3 repre- 
sents a very strong Pauli repulsion which increases with 
the fermion number Np ■ The purpose of this study is to 
show that a sufficiently strong attractive boson-fermion 
coupling term Np B \xB | 2 can overcome this repulsion and 
form the bright solitons. 



III. NUMERICAL RESULT 

We solve the coupled mean-field-hydrodynamic equa- 
tions (|2.15() and l|2.16|l for bright solitons numerically 
using a time-iteration method based on the Crank- 
Nicholson discretization scheme elaborated in Ref. [2l| . 
We discretize the mean-field-hydrodynamic equation us- 
ing time step 0.0005 and space step 0.025. 

We performed the time evolution of the set of equa- 
tions (|2.15() and (|2.16() introducing harmonic oscillator 
potentials y 2 in these equations and setting the nonlin- 
ear terms to zero: N BB = N B p — Np B — Nff = and 
starting with the eigenfunctions of the linear harmonic 
oscillator problem, e.g., with x_b(j/,t) = XF(y,i~) — 
7T" 1 / 4 exp(— y 2 /2) exp(— ir). The introduction of the ex- 
tra harmonic oscillator potential in these equations only 
aids in starting the time evolution with an exact analytic 
solution. In the end the harmonic oscillator potentials 
will be set equal to zero and will have no effect on the fi- 
nal wave function for solitons. During the course of time 
evolution the nonlinear terms are switched on very slowly 
and the resultant solutions iterated until convergence was 
obtained. Then the time evolution is continued and the 
harmonic oscillator potential terms in both bosonic and 
fcrmionic equations are slowly switched off by reducing 
the y 2 term to zero in 10000 steps of time evolution. Then 
the resultant solutions are iterated about 50000 times for 
convergence without any harmonic oscillator potential. If 
converged solutions are obtained, they correspond to the 
required axially-free bright solutions in the absence of 
any axial potential. In our numerical investigation as in 
the theoretical study of Refs. 0,0] we use u> = 2tt x 100 
Hz, and take m B as the mass of 87 Rb. Consequently, the 
unit of length I w 1 /im and unit of time 2/lo sa 3 ms. 

First we solve Eqs. (|2~l"5l) and [ETTtJjl with N F = N B = 
1000, a BB = 5 nm and a B p — —20 nm. This value of 
a B p is the experimental scattering length of Rb atoms 
|l3| . and a B F = —20 nm is the experimental scattering 
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FIG. 1: (Color online) The stationary function |</>i(2,i)| for 
axially-free bosonic (dotted line) and fermionic (full line) 
bright solitons vs. z for N F = N B = 1000, a B B = 5 nm, 
aBF ~ —20 nm, harmonic oscillator length / w 1 /j,m and 
v = 0. The arrows in dotted and full lines indicate the 
bosonic and fermionic axes, respectively. The nonlinearity 
parameters are Nbb = 20, Nbf = —160, N F b ~ —160, and 
N FF = 274.6. 
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FIG. 2: (Color online) The stationary function |0i(z,t)| for 
axially-free bosonic (dotted line) and fermionic (solid line) 
bright solitons vs. z for N F = 1000, N B = 10000, a B B = -1 
nm, asF = —1.875 nm, harmonic oscillator length I ~ 1 /im 
and v = 0. The arrows in dotted and full lines indicate the 
bosonic and fermionic axes, respectively. The nonlinearity 
parameters are Nbb = —40, N FF ~ —15, N FF ~ —150, and 
N FF = 274.6. 



length of the Rb-K system [|| . With these parameters 
the nonlinearities in Eqs. (|2.15l) and H2.16fl are N BB = 
20, N BF = -160, N FB = -160, and Nff = 274.6. 
The converged bright solitons are plotted in Fig. 1. In 
this case the fermionic and bosonic functions, 4>f and 
<f> B , respectively, have similar spatial extentions. It is 
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FIG. 3: (Color online) The stationary function |0i(.z,t)| for 
axially-free bosonic (dotted line) and fermionic (solid line) 
bright solitons vs. z for iV F = 1000, N B = 10000, a BB = 0.5 
nm, clef = —3.75 nm, v = and harmonic oscillator length 
I ss 1 fim. The arrows in dotted and full lines indicate the 
bosonic and fermionic axes, respectively. The nonlinearity 
parameters are Nbb = 20, Nbf = —30, Nfb = —300, and 
N FF = 274.6. 



possible to have solitons with different extensions in space 
by varying the parameters of the system. We took the 
experimental values for the scattering lengths in Fig. 1. 
However, the scattering length can be manipulated in the 
boson-fermion 6 Li- 23 Na and 40 K- 87 Rb systems near the 
recently discovered Feshbach resonances in them by 
varying a background magnetic field. Thus by varying 
the scattering length and the number of atoms we could 
arrive at different values of nonlinearity than in Fig. 1. 

To simulate a different situation of nonlinearity pa- 
rameters next we take a BB = — 1 nm, a B p = —1.875 
nm, N F = 1000, and N B = 10000. so that N BB = -40, 
N BF = -15, N FB = -150, and Np p — 274.6. In this 
case the profiles of the bright solitons shown in Fig. 2 are 
very different from those in Fig. 1. In Fig. 1 both the 
solitons are localized to a small region in space, whereas 
in Fig. 2 only the bosonic soliton is localized to a small 
region in space whereas the fermionic soliton extends to 
a very large region of space. 

Next we consider N F = 1000, N B = 10000, a BB = 
0.5 nm and a B F = —3.75 nm, so that the nonlinearity 
parameters are N BB = 20, N B p — —30, Np B = —300, 
and Npp = 274.6. The profiles of the solitons in Fig. 3 
are different from those in Figs. 1 and 2. In this case the 
bosonic function extends over a longer region in space 
than the fermionic function. In Figs. 1 and 2 it was the 
fermionic function that extends over a longer region in 
space. 

In Fig. 2 the nonlinearity Np B appearing in the 
fermion equation is less attractive compared to that in 
Figs. 1 and 3. Hence the resultant nonlinear interaction 
in the fermion component is less attractive and hence the 



|<yz,t)| (ixm- 1/2 ) 




z(nm) 



30 
20 
10 

(b) 



|(j) F (z,t)| (um 1/2 )J 



Fermion 




zQxm) 



FIG. 4: (Color online) The function \<f>i(z, t)\ for (a) bosonic 
and (b) fermionic bright solitons vs. z and t for the solitons 
of Fig. 3. At t = N F = 1000, N B = 10000, a BB = 0.5 nm, 
aBF = —3.75 nm, v — and harmonic oscillator length / ~ 1 
fim. The nonlinearity parameters at t = are Nbb ~ 20, 
N BF = -30, N FB = -300, and N FF = 274.6. At t = 100 
ms (marked by arrows) the bright solitons are set into small 
breathing oscillations by suddenly jumping the nonlinearities 
Nbf and Nfb to —33 and —330, respectively. 



fermionic soliton extends to a large distance in space. If 
the attraction in Np B is further reduced the fermionic 
soliton ceases to bind. 

Hence by manipulating the parameters one could have 
different situations of localization of the solitons. One 
could cither have both the solitons extending up to sim- 
ilar distance in space as in Fig. 1 or one of the solitons 
extending over a longer region in space as in Figs. 2 and 
3. It is worth emphasizing that in the fermionic equation 
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FIG. 5: (Color online) The root-mean-square size {z) rms of 
the bosonic (dotted line) and fermionic (full line) solitons of 
Fig. 4 vs. time. 



(|2.1(il) the diagonal nonlinearity Npp is repulsive, hence 
the binding solely comes from the attractive off-diagonal 
nonlinearity Nfb corresponding to an attractive boson- 
fermion interaction. Hence for a fermionic soliton to ap- 
pear Nfb and asF are always taken to be negative or 
attractive. Consequently, Nbf is also negative. Finally, 
terms Nbb and cibb could be either positive or negative. 
When these terms are positive or repulsive the bosonic 
solitons are formed due to an attractive or negative Nbf, 
as in Fig. 3. 

Next we study the stability of the bright solitons nu- 
merically. We consider the soliton of Fig. 3 and dur- 
ing time evolution we suddenly jump at t = 100 ms the 
nonlinearity Nbf from —30 to —33 and the nonlinear- 
ity Npb from —300 to —330. This can be achieved by 
manipulating a background magnetic field near a Fesh- 
bach resonance |22j in the boson-fermion interaction and 
thus varying the boson-fermion scattering length by 10%. 
Due to the sudden change in the nonlinearity the bosonic 
and fermionic bright solitons are set into stable breath- 
ing oscillation. The evolution of the wave function profile 
of the two solitons are shown in Figs. 4. The solitons 
are found to execute stable non-periodic breathing os- 
cillation. In Fig. 5 we plot the root-mean-square size 
{z)rms of the bosonic and fermionic solitons of Fig. 4 as 
a function of time. The breathing oscillation of the two 
solitons after the perturbation is applied results in the 
stable non-periodic oscillation of the root-mean-square 
sizes illustrated in Fig. 5. As, after applying the pertur- 
bation, the boson-fermion attraction has been increased 
this corresponds to a reduction in the root-mean-square 
size (z) rms as we find in Fig. 5. The steady propagation 
of the solitons in Figs. 4 and the stable oscillation of their 
root-mean-square sizes in Fig. 5 after the perturbation 
is applied demonstrate the stability of the solitons. 

Finally, we consider the fermionic bright solitons 



formed on a periodic optical-lattice potential. For that 
purpose we include in Eqs. i|2.15[l and (|2.1rj|) the follow- 
ing optical-lattice potential formed by a standing-wave 
laser beam [l7| 



Vqp = sin 2 (2^/ A), 



(3.1) 



where Vq is the strength, and A is the wave length of the 
laser. In our calculation we take Vb = 100 and A = ir/2. 
To solve Eqs. I|2.15|l and i|2.16[l with this optical-lattice 
potential and desired nonlinearities, we again start the 
time evolution with the solution of the linear oscillator 
problem. In the course of time evolution the nonlinear- 
ities and the optical-lattice potential are slowly intro- 
duced and eventually the harmonic oscillator potential 
is slowly removed. Then the final solutions are iterated 
for convergence. The resultant soliton wave functions are 
plotted in Fig. 6 for N F = 1000, N B = 10000, a BB = 0.3 
nm, asF — —2.375 nm. The optical-lattice potential in- 
troduces modulations in the solitonic wave function. For 
the parameters of Fig. 6 the modulations are more promi- 
nent on the fermionic soliton than the bosonic one. By 
changing the parameters it is possible to have modula- 
tions on the bosonic soliton as well. 
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FIG. 6: (Color online) The stationary function |<^i(z,t)| for 
bosonic (dotted line) and fermionic (full line) bright soli- 
tons vs. 2 for N F = 1000, N B = 10000, a B B = 0.3 nm, 
aBF = —2.375 nm, v = and harmonic oscillator length 
I ~ 1 fim in the presence of the optical-lattice potential 
V(y) — 100 sin 2 (4j/). The arrows in dotted and full lines indi- 
cate the bosonic and fermionic axes, respectively. The nonlin- 
earity parameters are Nbb = 12, Nbf = —19, Nfb = —190, 
and Nff — 274.6. 



IV. SUMMARY 

We use a coupled set of time-dependent mean-field- 
hydrodynamic equations for a boson-fermion mixture to 
study the formation of fermionic bright soliton in a DFG 
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as a stationary state. In this study we take the boson- 
boson interaction to be both attractive and repulsive and 
the boson-fermion interaction to be attractive. An at- 
tractive boson-fermion interaction is necessary for the 
formation of a fermionic bright soliton as the diagonal 
nonlinearity Nff in the fermion-fermion system is al- 
ways repulsive. 

In the present study we demonstrate that stable soli- 
tons can be formed in coupled NLS equations for the 
boson-fermion mixture with repulsive diagonal nonlin- 
earities and attractive off-diagonal nonlinearities above 
some cut-off values. In another study [2^ we showed 
the possibility of the formation of bright solitons in cou- 
pled bosonic condensates with intraspecies repulsion sup- 
ported by interspecies attraction. The stability of the 
present fermionic and bosonic solitons is demonstrated 
through their sustained breathing oscillation initiated by 
a sudden jump in the boson-fermion scattering length. 
Bright solitons have been created experimentally in at- 
tractive BECs in three dimensions in the presence of ra- 



dial trapping only without any axial trapping [l5j|. I n 
view of this fermionic bright solitons can be observed in 
the laboratory in the presence of radial trapping only in 
a mixture of a DFG and BEC. We also suggest the pos- 
sibility of the formation of fermionic solitons on a peri- 
odic optical-lattice potential. In the present investigation 
we used a set of mean-field equations for the DFG-BEC 
mixture. A proper treatment of the DFG should be per- 
formed using a fully antisymmetrized many-body Slater 
determinant wave function. However, we do not believe 
that the present conclusion about the existence of robust 
fermionic solitons in a DFG-BEC mixture to be so pecu- 
liar as to have no general validity. 
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